#ifndef MATH_OP_H
#define MATH_OP_H

#define EPSINON (0.0000001)

bool is_close(const double &f1, const double &f2, const double &epsinon=EPSINON)
{
    const double diff = f1 - f2;
    return (diff >= -epsinon) && (diff <= epsinon);
}


#endif // MATH_OP_H
